Cosmic Microwave Background and Supernova Constraints on Quintessence: 
Concordance Regions and Target Models 
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We perform a detailed comparison of the Wilkinson Microwave Anisotropy Probe (WMAP) mea- 
surements of the cosmic microwave background (CMB) temperature and polarization anisotropy 
with the predictions of quintessence cosmological models of dark energy. We consider a wide range of 
quintessence models, including: a constant equation-of-state; a simply-parametrized, time-evolving 
equation-of-state; a class of models of early quintessence; scalar fields with an inverse-power law po- 
tential. We also provide a joint fit to the CBI and ACBAR CMB data, and the type la supernovae. 
Using these select constraints we identify viable, target models for further analysis. 

PACS numbers: 98.80.-k 



The precision measurement of the cosmic mi- 
crowave background (CMB) by the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite P, Q represents a 
milestone in experimental cosmology. Designed for pre- 
cision measurement of the CMB anisotropy on angular 
scales from the full sky down to several arc minutes, this 
ongoing mission has already provided a sharp record of 
the conditions in the Universe from the epoch of last 
scattering to the present. In light of this powerful data 
0, 0, IH B , we must consider anew our cosmological 
theories. 

We aim to use the WMAP results to test cosmological 
theories of the accelerating Universe — to seek clues to 
the nature of the dark energy. Despite the absence of a di- 
rect dark-energy interaction with our baryonic world, the 
CMB photons provide a probe of the presence of the dark 
energy, complementary to the type la supernovae. Via 
the integrated Sachs- Wolfe effect on large angular scales, 
the geometric optics of the last-scattering sound horizon 
on degree scales, and the pattern of acoustic oscillations 
on smaller angular scales, we expect the CMB to reveal 
information about the dark energy density, equation-of- 
state, and behavior of fluctuations. These effects are il- 
lustrated in Figure n 

In this article we test the quintessence hypothesis, that 
a dynamical, time-evolving, negative pressure, inhomo- 
geneous form of energy dominates the cosmic energy 
density and is responsible for the cosmic acceleration 
milllllllllllllll- To be precise, we carry out 
an extensive analysis of the CMB anisotropy and mass 
fluctuation spectra for a wide range of quintessence mod- 
els. These models are: (Ql) models with a constant 
equation-of-state, w, including w < —1; (Q2) models 
with a simply-parametrized, time-evolving w; (Q3) early 
quintessence models, with a non-negligible energy density 
during the recombination era; (Q4) trackers described by 
a scalar field evolving under an inverse-power law poten- 
tial. 

The suite of parameters describing the cosmologi- 
cal models are split into spacetime plus "matter sec- 
tor" variables, 0m, and separate quintessence parame- 
ters, 9q. The spacetime and "matter sector" of the 
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FIG. 1: The pattern of CMB anisotropy can reveal infor- 
mation about the quintessence abundance (Oq), equation-of- 
state (w), and behavior of fiuctuations (S). The three curves 
are examples of constant equation-of-state models which dif- 
fer little by eye, but are distinguished by the data. The red 
{w = —0.5) and blue ( — 1.2) curves are both low-x^ CMB- 
indistinguishable, but distinct with respect to SNe. The black 
curve (—0.8), although it is consistent with the SNe data and 
matches the location and height of the first acoustic peak 
determined by WMAP P, is rejected by the CMB at the 
3(T-level. 



quintessence models are specified by the parameter set 
6m = {^bh'^, ^cdmh"^, h, Us, As, r^}. In order, these 
are the baryon density, cold dark matter density, bub- 
ble parameter, scalar perturbation spectral index, scalar 
perturbation amplitude, and optical depth. In this in- 
vestigation we restrict our attention to spatially-flat, cold 
dark matter models with a primordial spectrum of nearly 
scale-invariant density perturbations generated by infla- 
tion. 

The quintessence parameters vary from model to 
model. For the simplest family of models, with a constant 
equation-of-state, we need only to specify 9q — {w}. For 
models which feature more realistic time-evolution of the 
quintessence, which may include a non-negligible frac- 
tion of quintessence at early times, more parameters are 
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required, e.g. 9q = {w, dw/da, ...}, to characterize the 
impact on the cosmology in general and the CMB in par- 
ticular. 

Our analysis method is as follows: (i) compute the 
CMB and fluctuation power spectra for a given cosmo- 
logical model; (ii) compute the relative likelihood of the 
model with respect to the experimental data; (iii) assem- 
ble the likelihood function in parameter space to deter- 
mine the range of viable quintessence models. For step 
(i) we use both a version of CMBfast niodified for 
quintessence, as well as the newly-available CMBEASY 
[l^. For step (ii) we supplement the WMAP data with 
the complementary ACBAR ^ and CBI-MOSAIC data 
[Tsl IT^ (using the same bins as Ref. IE H), in addi- 
tion to the current type la SNe data IM l2ll I2I E^ . 
Certain other constraints, such as the HST Key Project 
measurement of Hq \2m or the limit from Big Bang Nu- 
cleosynthesis on rifc/i^ [23 through the deuterium abun- 
dance measurement are satisfied as cross-checks. It is 
remarkable that such agreement can be found between 
such diverse phenomena. Our focus in the following in- 
vestigation, however, is primarily on the CMB and SNe. 
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FIG. 2: The one-parameter family of best-fit models, which 
exploit the geometric degeneracy of the CMB anisotropy pat- 
tern, is shown as the thick, red curve in the w — h plane. We 
have explored models in a six-dimensional cylinder in the pa- 
rameter space surrounding this "best-fit line." The HST-Key 
Project la measurement of the Hubble constant is shown by 
the shaded band. 

Ql: We have analyzed the cosmological constraints on 
the simplest model of quintessence, characterized by a 
constant equation-of-state, w. We have used the equiv- 
alence between a scalar field if with potential V{ip) and 
the equation-of-state w in order to self-consistently eval- 
uate the quintessence fluctuations. For the range w < —I 
we employ a k essence model, keeping the sound speed 
(actually, this is duj'^/dk'^) fixed at = 1. Since this 
model introduces only one additional parameter beyond 
the basic set of spacetime plus matter sector variables, 
we adopt a simplistic grid-based search for viable mod- 
els. The acceptance criteria for the Ql models is based 
on a Ax^-test. The results of our survey of Ql models 
are shown in Figures [3121 We have exploited the de- 
generacy of the CMB anisotropy pattern among models 
with the same apparent angular size of the last scat- 
tering horizon ^26]. Hence, there is a family of models 
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FIG. 3: The constraints on constant equation-of-state models 
due to CMB (WMAP, ACBAR, CBI) and type la supernovae 
(Hi-Z, SCP) are shown. The starting point for our parameter- 
search, the family of CMB-degenerate models, is shown by the 
thick, black line. 

with ribh"^ = 0.023, ncdmh"^ = 0.126, = 0.97, and 
characterized by pairs {w, h} having (nearly) indistin- 
guishable CMB anisotropy patterns. The pairs {w, h} 
are shown in Figure |21 and all represent quintessence 
models with = 1429 for the WMAP temperature- 
temperature and temperature-polarization data — a one- 
parameter family of best-fit models. From this starting 
point, we explored over 6 x 10^ models distributed on 
a grid filling a 6-cylinder around the best-fit line, vary- 
ing {rifjft,^, Q.cdm.h'^ , h, rig, Tr\ at intervals in w. For each 
model we evaluate the likelihood relative to WMAP, as 
well as the complementary ACBAR and CBI-MOSAIC 
data. The 2a boundary, based on a Ax^ test for six 
degrees of freedom is shown in Figure |31 We have also 
evaluated the constraint in the w — plane for the 
combined Hi-Z/SCP type la supernova data set, show- 
ing the 2(T region based on a Ax^ test for two degrees 
of freedom. Our basic conclusion from the overlapping 
constraint regions is that there exist concordant models 
with -1.25 < w < -0.8 and 0.25 < < 0.4. We 
have identified four sample models in Tabled for further 
analysis. 



model 


w 


h 




Ql.l 


-0.82 


0.630 


0.84 


Q1.2 


-1.00 


0.682 


0.89 


Q1.3 


-1.18 


0.737 


0.96 


Q1.4 


-1.25 


0.759 


0.97 



TABLE L Sample best-fit models with a constant equation-of- 
state (Ql). All models have n^h'^ = 0.023, n^dmh^ = 0.126, 
Hs — 0.97, and Tr = 0.11. 

Our search of the parameter space for the remaining 
models is based on a Bayesian approach, using a Monte- 
Carlo Markov-Chain (MCMC) search algorithm to iden- 
tify the best cosmological models. The end product is a 
realization of the posterior probability distribution func- 
tion on the parameter space p7j . Our approach is similar 
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(a) Q2: monotonic evolution w{a) 



(b) Q3: leaping kinetic quintessence 
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(e) Q3: leaping kinetic quintessence 



S 

a 




0.05 



(f) Q4: inverse-power law 



FIG. 4: The results of our MCMC search of the multi-dimensional parameter search, for models Q2-4, are illustrated in the 
figures above. In all cases, we have marginalized over the suppressed parameters. The solid lines indicate the 1, 2, 3cr contours 
based on comparison with the CMB (WMAP, ACBAR, CBI) and type la supernovae (Hi-Z, SCP). 



to the procedure described in 0, whereby the MCMC 
makes a "smart" walk through the parameter space, ac- 
cepting or rejecting sampled points based on a running 
criteria. For each of Q2-4, after some experimentation we 
found it practical to use four independent chains in or- 
der to monitor convergence and mixing according to the 
criteria of Ref. Each such chain explored 3 x IC* 
models. We found it advantageous to use a dynamical 
stepsize factor r S [0.1,10]. In each step, all parame- 
ters are varied according to a Gaussian with width rai, 
where cr^ is an initial stepsize for the random walk in the 
direction of the i-th parameter. The stepsize factor r is 
adjusted such that the chain neither rejects nor accepts 
too frequently. Certain parameters, such as wq and fig in 
the early quintessence models, have physical boundaries 
and it is therefore important that our MCMC does not 
lead to a bias in the sampled distribution close to the 
boundaries. As the likelihood in the MCMC approach 
is given by the frequency with which the chain passes 
through a given region of parameter space, we would be 
in danger of depleting the region close to the boundary 



if we were to simply discard any steps which cross the 
boundary. That is, the likelihood in the parameter region 
close to the boundary would be artificially suppressed. 
If such a (virtual) jump in a parameter direction i oc- 
curs, we use the following strategy. We try making a 
gaussian-weighted step but only towards the boundary. 
If this fails after ten attempts, then we start over again 
by making a gaussian step in any direction. By trying 
ten times, however, it is very likely that we do make the 
step successfully, if we are at least one ai away from the 
boundary. If, however, we are very close to the boundary 
already and our step would have carried us deep into the 
unphysical region, then it is sensible to make the step 
away from the current point, for the chain didn't intend 
to stay close to the boundary anyhow. To give us confi- 
dence in the reliability of this algorithm, we have checked 
by synthetic distributions that this strategy is very well 
suited for sampling the likelihood, even in cases where a 
physical boundary cuts off the parameter space. 

Q2: We have examined quintessence models with an 
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equation-of-state that evolves monotonically with the 
scale factor, as w{a) = Wq + {I — a)wi. For this case, the 
parameters are simply = {wq, wi}. This parametriza- 
tion has been shown to be versatile in describing the late- 
time quintessence evolution for a wide class of scalar field 
models [2^ . Based on the degeneracy of models found for 
Ql, we expect to find a two-dimensional family of equiv- 
alent best-fit models with the same apparent angular size 
of the last scattering horizon, occupying a plane in the 
{wq, Wl, h} space. There are three ways in which this 
plane is pared down: Firstly we confine w > —1 at all 
times. Secondly, the transition from wq to wq + wi takes 
place at low redshift, z < 1, so that high redshift su- 
pernovae restrict wq, wi for these models. Thirdly, this 
parametrization allows for models in which the dark en- 
ergy is non-negligible at early times, which influences the 
small-scale fluctuation spectrum. The first two consider- 
ations yield wq < —0.75 at the 2(7 level, marginalizing 
over the suppressed five dimensional parameter space, as 
illustrated in Figure 4(a) There, the shapes of the con- 
tours indicate that current data can only distinguish be- 
tween fast {dw/da > 0.5) and slow evolution of w(a), and 
offer only a weak bound on wi. However, in terms of Jig, 
our third consideration gives a tight upper bound on the 
quintessence density during recombination. As shown 
in Figure |4(d]] < 0.03 at the 2cr level. Three tar- 



get models, with significantly different equation-of-state 
evolution dw/da, are given in Table UTI for future investi- 
gations. 



model 


wo 


Wl 


h 




crs 


Q2.1 


-0.93 


0.43 


0.66 


4 X 10"'= 


0.77 


Q2.2 


-0.99 


0.68 


0.64 


1 X 10~* 


0.78 


Q2.3 


-0.92 


0.62 


0.62 


1 X 10~* 


0.73 



more general parameterization, allowing for independent 
Wo and steepness of transition can be found in js^l .) Since 
fig is not tied as closely to the expansion rate sampled by 
the supernovae, compared to case Q2, the result is the 
weaker constraint fig < 0.1, as shown in Figure 4(b) 
Although the limit of a cosmological constant can be ap- 
proached in this model, the presence of early quintessence 
will then require a sharp transition in the equation-of- 
state in order to reach w —1. In addition to the fact 
that such models lose the early tracking behavior and 
instead require some degree of fine-tuning, there is the 
practical consideration that the sharp transition leads to 
some numerical instability in our code. To avoid this 
problem, we restrict w > —0.97 as can be seen in Fig- 
ure 4(b) We believe that we have succeeded at imple- 
menting an MCMC algorithm that does not artificially 
distort the probability distribution near the parameter- 
space boundaries, based on our experimentation with 
synthetic distributions with boundaries. Next, because 
early quintessence suppresses the growth of fluctuations 
on small scales compared to large scales, we find that 
comparable fluctuation spectra can be achieved by mak- 
ing a trade-off between Us and fig. As shown in Fig- 
ure |4(e)| slight suppression of small-scale power can be 
accomplished either by a tilt towards the red, < 1, 
or a rise in fig . Since the effect of early quintessence on 
the small-scale fiuctuation power spectrum closely mim- 
ics a running spectral index, we have not introduced 
drts/dln/c as an additional parameter, which would be 
highly degenerate with fig jSlj]. We expect that im- 
proved measurements of the second and third acoustic 
peaks will tighten the constraint on Slg and sharpen the 
degeneracy in the rig — f2g plane. Target models, with 
significantly different values of early quintessence abun- 
dance, are given in Table UTTl 



TABLE IL Sample best-fit models with a monotonically- 
evolving equation-of-state (Q2). Although a range of param- 
eters give equivalently good fits to the observational data, we 
have selected this sample with {ftbh^ , ftcdmh^, ris, r^} — 
{0.023, 0.11, 0.98, 0.16} (Q2.1), {0.023, 0.11, 0.98, 0.15} 
(Q2.2), {0.023, 0.12, 0.98, 0.08} (Q2.3), Entries for Q}^ and 
as, are the resulting values based on the other parameters. 



Q3: We have examined models of leaping kinetic 
quintessence, a scalar field evolving under an exponen- 
tial potential with a non-canonical kinetic term that un- 
dergoes a sharp transition at late times, leading to the 
current accelerated expansion |33|- At early times the 
field closely tracks the cosmological background with 
w = during matter domination, appearing as early 
quintessence [31] before undergoing a steep transition 
towards a strongly negative equation-of-state by the 
present day. The steepness of the transition in w for 
a leaping kinetic model is directly connected to the 
equation-of-state wq today. Such models can therefore be 
characterized by the parameters 6q — {{1^q,wq}, where 
rig is the quintessence density during recombination. (A 



model 


wo 


nils 


h 


A 


Wis 


0-8 


Q3.1 


-0.94 


0.006 


0.69 


0.0026 


-0.27 


0.89 


Q3.2 


-0.91 


0.024 


0.70 


-0.0028 


-0.21 


0.81 


Q3.3 


-0.93 


0.043 


0.71 


-0.0070 


-0.19 


0.85 



TABLE IIL Sample best-fit leaping-kinetic quintessence 
models (Q3). The parameters A is an equivalent descrip- 
tion of fig, which along with the averaged w at last scat- 
tering, wj^, can be used more easily with equations 2-4 of 
Ref. pHl to generate the time-evolution of these models. The 



non-quintessence parameters are {Q.bh , ilcdmh 



{0.022, 0.112, 0.96, 0.12} (Q3.1), {0.023, 0.116, 1.0, 0.16} 
(Q3.2), {0.024, 0.119, 1.04, 0.26} (Q3.3). 



Q4-: Finally, we have examined tracker models of 
quintessence. Inverse-power law (IPL) models are the 
archetype quintessence models with tracking property 
and acceleration 0, IH, Q . The potential is given by 
V cx ip~" , where the constant of proportionality is deter- 
mined by Og. In certain supersymmetric QCD realiza- 
tions of the IPL I33 , a is related to the numbers of color 
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and flavors, and can take on a continuous range of values 
a > 0. For a ^ 0, however, inverse-power law mod- 
els behave more and more like a cosmological constant. 
Using earlier data, a has been constrained to a < 1.4 
[sgLIs^ . although h — 0.65 has been fixed in those anal- 
yses. Keeping h free, a more conservative value of a < 2 
[38| was inferred. From our analysis, we see that the 2a 
bound has not changed dramatically. The 2(T-bound with 
a < 1 — 2 is consistent with values of h within the range 
determined by the HST, as seen in Figure [4(c)| In Fig- 
ure 4(f) we plot the likelihood contours in the ilmh^ — ct 
plane: our results agree with the best fit at ^Imh^ = 0.149 
for a = Oorw = — 1, but show a tolerance for a wider 
range for < a < 2. That is, the additional degree of 
freedom in a means that the matter density for the IPL 
model is not as well-determined from the peak position 
as compared to the A model. However, to maintain 
the peak ai £ = 220, we observe that ri„i/i^ decreases 
slightly as a increases. We might have inferred the re- 
sults for the IPL based on the constant equation-of-state 
models: pairs of {a, h} can equivalently determine a 
family of models with degenerate CMB anisotropy pat- 
terns, since the differences in the late-ISW for this model 
compared to Ql make only a small contribution to the 
overall x^- Furthermore, since IPL quintessence can be 
modeled by the appropriate choice of the Q2 parameters, 
then improved sensitivity to dw / da is required to tighten 
the constraints here. While a is tightly constrained, our 
take-away is that IPL models with 0.25 < S 0-4 
remain viable. Target models are given in Table Hvl 



model 


a 


h 


0-8 


Q4.1 


0.1 


0.68 


0.90 


Q4.2 


0.2 


0.68 


0.85 


Q4.3 


0.8 


0.68 


0.82 



TABLE IV: 
els (Q4). 



Sample best-fit IPL quintessence mod- 
The non-quintessence parameters are 
Us, Tr] = {0.023, 0.122, 0.97, 0.13} (Q4.1), 



{0.023, 0.116, 0.97, 0.14} (Q4.2), {0.024, 0.102, 1.0, 0.23} 
(Q4.3). 

This work provides a capsule summary of the vi- 
able quintessence dark energy models, based on two of 



the tightest constraint methods, using the CMB and 
SNe. Our study advances beyond past investigations 
EllliliSlilEllillillil by treating a wide class 
of quintessence models with the powerful weight of the 
WMAP data. We have considered four classes of models 
which cover the most basic quintessence scenarios, in- 
cluding a versatile parametrization, as well as the best 
motivated and most realistic scenarios based on our cur- 
rent understanding of particle physics. Absent from our 
survey are k essence models, and dark energy models with 
a coupling to other matter fields or gravity. In the former 
case, since the sound speed of fluctuations varies with 
time in these models, however, we can make a simple 
distinction with quintessence models with an underlying 
scalar field, wherein the propagation speed is equal to the 
speed of light. (See Refs. ^St 43 for analysis of these 
models with respect to CMB anisotropy.) For the latter 
case we refer to Ref. plj for specific coupled models. We 
also note that the mass fluctuation power spectrum is an 
important cosmological constraint which we have omit- 
ted at this stage, primarily because it constrains energy 
density rather than pressure (although there are excep- 
tions), a chief feature distinguishing dark energy from 
dark matter. (The constraint curves obtained in Ref. 
|4^ . e.g. Figure 3, are consistent with, but do not de- 
cisively pare down the parameter regions of our current 
results.) Furthermore, the analysis of mass power spec- 
trum observations will have to take into account possible 
influences of a time-dependent w and early quintessence, 
which we put off for later investigation. Overall, we have 
simulated in detail more than 400,000 individual cosmo- 
logical models. The stored spectra and parameter-space 
likelihood functions will be used to evaluate additional 
constraints that can offer clues to the behavior of the 
dark energy. The set of thirteen target models, listed 
in Tables IIIIVI which includes a revised concordance A 
model 50] , should be of use by researchers to further test 
the quintessence hypothesis. 
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